Gate driven adiabatic quantum pumping in graphene 
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We propose a new type of quantum pump made out of graphene, adiabatically driven by oscil- 
lating voltages applied to two back gates. From a practical point of view, graphene-based quantum 
pumps present advantages as compared to normal pumps, like enhanced robustness against thermal 
effects and a wider adiabatic range in driving frequency. From a fundamental point of view, apart 
from conventional pumping through propagating modes, graphene pumps can tap into evanescent 
modes, which penetrate deeply into the device as a consequence of chirality. At the Dirac point the 
evanescent modes dominate pumping and give rise to a universal response under weak driving for 
short and wide pumps, even though the charge per unit cycle in not quantized. 



I. INTRODUCTION 

Externally driven quantum devices provide an attrac- 
tive setting to transfer charges between electronic reser- 
voirs [IHI]. Conventional realizations of such quan- 
tum pumps are based on patterned semiconductor het- 
erostructures, and commonly utilize the adiabatic or non- 
adiabatic loading and unloading of electrons into well- 
defined individual quantum states, where energy gaps are 
provided by Coulomb interactions or by size quantization 

ishu]. 

Graphene, consisting of an atomically thin sheet 
of carbon atoms [13 , provides a novel platform for 
two-dimensional electronic transport which combines a 
regime of nominally vanishing charge carrier density (at 
the Dirac point) with radically different confinement 
properties. These properties are intimately linked to 
the chiral nature of the charge carriers, which suppresses 
backscattering at interfaces and results in the so-called 
Klein paradox by which charge carriers are difficult to 
confine [T4HT6] . Here we show that chirality and vanish- 
ing carrier density conspire to yield a universal dimen- 
sionless pumping efficiency at the Dirac point, which re- 
markably corresponds to a non-quantized charge per unit 
cycle [17]. The pumping then is mediated by evanescent 
modes, which previously have been studied in the con- 
text of stationary transport [18-20 . We also analyze the 
behavior of graphene pumps away from the Dirac point 
jTTJ [2TJ [22], which is dominated by propagating modes, 
and compare it to that of normal quantum pumps. In 
the second part of this paper we turn to issues of prac- 
tical implementations. We find that graphene's spectral 
properties enhance the response to external driving by 
gates, increase the robustness against thermal effects, 
and result in a wider adiabatic range in driving frequency. 
Graphene pumps therefore offer a number of potential 
advantages which may prove useful for practical applica- 
tions. 

This work is organized as follows. Section [IT] intro- 
duces a model of a quantum pump driven by two electro- 
static gates, which we analyze in the scattering approach 
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FIG. 1. (Color online) Top panel: Graphene quantum pump 
where two gates at voltages Vi(t) and Vb(£) induce a peri- 
odic time-dependent onsite energies U\(t) and U2(i) in the 
two halves of a graphene flake. We investigate the resulting 
charge transport between two contact electrodes, separated 
by a distance 2L. Bottom panels: Instantaneous onsite po- 
tential at a point during the pumping cycle (solid thick line), 
including the potential in the contact region for highly doped 
contacts (left, Ul,r — ► —00) and undoped contacts (right). 
The dashed line shows the potential at the working point (in 
the absence of driving) . 



to adiabatic quantum pumping. Results for normal and 
graphene-based pumps are presented in Sec. |III[ where we 
formulate them in terms of the onsite potential driving. 



In Sec. [TV] we apply these results to the experimentally 
relevant case of driving via gate voltages. Section [V] con- 
tains estimates of the maximal pumping performance of 
realistic graphene devices and their normal counterparts. 
We conclude in Sec. [VI] with a summary of the practical 
advantages of graphene-based quantum pumps. 



II. 



MODEL AND FRAMEWORK 



In order to pump an average current between two reser- 
voirs that are kept at the same bias, it is necessary to vary 
the scattering properties of the pump region (the system 
between the reservoirs) periodically over time. This is 
achieved by driving the system with some sort of external 
forces that modify the parameters of the system, hence 
the name 'parametric pumping'. The driving forces can 
be induced by any external parameter (like a magnetic 
field, the transparency of the contact between the pump 
and the reservoir, the potential energy, etc.), as long as 
their variation modifies the scattering matrix of the sys- 
tem. When the temporal variation of these parameters 
is sufficiently slow (as compared to the dwell time of the 
carriers within the pump region), the pump operates in 
the adiabatic regime. To each traversing electron, the 
system then is approximately static, and the scattering 
amplitudes become insensitive to the explicit time de- 
pendence of the driving. In this case, the pumped charge 
increases proportionally to the frequency of the external 
driving, and at least two parameters of the system must 
be driven out of phase. The direction of the pumped cur- 
rent depends on the specific way the scattering phases de- 
pend on the varying parameters, and changes sign when 
the driving cycle is reversed. 

In this work we consider a new type of adiabatic quan- 
tum pump, made out of graphene, and compare it to the 
equivalent normal pump. We propose to externally drive 
the pump by the capacitative coupling to two adjacent 
gates. Such devices have already been experimentally in- 
vestigated for stationary transport, both for graphene as 
well as for 'normal' conductors (such as semiconducting 
quantum wells). Oscillating gate voltages V\(t) and V^(t) 
induce an oscillating shift U\{t) and [^(t) in the onsite 
potential of each of the two adjacent regions coupled to 
the gates. These two oscillating potential barriers have a 
width W and a length L, with a total pumping region of 
length 2L [23 . A sketch of the described setup, attached 
to two metallic contacts (left and right reservoirs), can 
be seen in the upper panel of Fig. [I] 

We characterize the unique features of quantum pump- 
ing in graphene by comparison of four different variations 
of the setup above: graphene pumps with either heavily 
doped contacts (bottom left panel of Fig. pi) or undoped 
contacts (bottom right panel) are contrasted with pumps 
where the graphene flake is replaced by a normal sys- 
tem accommodating an ordinary two-dimensional elec- 
tron gas. In the absence of driving, i.e., at the 'working 
point' of the pump (dashed line in the bottom panels of 
Fig. [I]), systems with undoped contacts have a uniform 
carrier density, equal in both the leads and in the central 
pumping region. In contrast, the charge carrier density 
in heavily doped leads is very large; this can be realized 
by metallic electrodes. The charge carrier density in the 
system can be uniformly shifted by another back (or top) 
gate covering both the pumping region and the leads (not 
shown in Fig. [T]), which allows to tune the working point 



of the device. For heavily doped contacts, the back-gate- 
induced change of the charge carrier density in the leads 
can be neglected. The different types of contacts allow 
us to isolate the influence of evanescent modes, which are 
only induced by the highly doped contacts. 

We follow the scattering approach to quantum pump- 
ing [2j [3], which relates the adiabatic transfer of non- 
interacting charge carriers to the parametric evolution of 
the instantaneous scattering matrix S(t). For two inde- 
pendent pumping parameters £ = {£1,^2} (the minimal 
requirement for adiabatic pumping) and single channel 
reservoirs, the charge transferred across the scattering re- 
gion from the left reservoir reduces to an integral over the 
area A enclosed by the driving path in two-dimensional 
parameter space, 
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where e is the electron charge, and L, R indicate left and 
right contacts. The scattering matrix can be expressed in 
terms of the transmission probability through the system 
T and the scattering phase <p = a — j3 given by 
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where r (r') and t (t f ) are the ^-dependent reflection and 
transmission amplitudes for electrons arriving from the 
left (right) reservoir. For wide quantum pumps depicted 
in Fig. [T] the number of channels is large. However, 
because of the quasi one-dimensional design where the 
onsite potential U is independent of the transverse co- 
ordinate y, the different channels remain decoupled. In- 
dexing the channels by a quantum number n, the total 
pumped charge is therefore given by a sum Q = J^ n Q n . 

Into the longitudinal x-direction, we model the poten- 
tial profile by two abrupt steps of equal length L and as- 
sume that the two driving parameters £ = {Ui(t), U2(t)} 
have zero average, with maximum amplitudes 8U\ and 
(5/72, respectively. Undoped contacts have the same Fermi 
momentum kp as the pumping region, while heavily 
doped contacts have a much larger Fermi momentum, 
which can effectively be taken as infinite [24] . 

For this set-up, the transmission probabilities T n and 
phases <\> n can be computed by simple wave matching, 
where graphene regions are described by the Dirac equa- 
tion (with Dirac velocity vp) [14 , and normal regions 
by the Schrodinger equation for quasiparticles of effec- 
tive mass m*. The results depend on the characteristic 
energy scale for the longitudinal confinement of the car- 
riers, given as E^ = hvp/L in the graphene case, and 
E^ = h 2 /(2m*L 2 ) for the normal conductor. Depend- 
ing on the amplitude Sui of the dimensionless driving 
energies Ui = Ui/El, the driving can be characterized as 
weak (Sui <C 1) or strong (5ui ^> 1). For weak driving, 



the charge 

Qn « d 2 u Q n (0) J d Ul du 2 = A u d 2 u Q n (0) (2) 

pumped in channel n becomes proportional to the small 
area A u ~ 5ui5u2 enclosed in parameter space by the 
driving cycle, wherein d 2 Q n (ui, 112) can be approximated 
by a constant. For strong driving the integral in Eq. (fiT) 
has to be performed numerically. 



while at x = L\ 



III. 



RESULTS AND DISCUSSION 



Let us first consider the case of graphene with heavily 
doped contacts. There are four regions in the scatter- 
ing problem: the left contact (x < —L), the region of 
the first barrier with onsite potential u\ (-L < x < 0), 
the region of the second barrier with onsite potential u^ 
(0 < x < L), and the right contact (x > L). The trans- 
port problem across the pump can be solved by match- 
ing the propagating and evanescent modes to the left and 
right of the interfaces separating different regions, which 
we obtain from the equation H^/ = (E — Ui)^ . Here H 
is the Dirac Hamiltonian H = vpa • p, where p is the 
momentum operator relative to the Dirac point and o\ 
are the Pauli matrices. The Dirac Hamiltonian acts on a 
two-component spinor, \I/ = (0a, </>b) T , representing the 
amplitude of the wavefunction of energy E on the two in- 
equivalent triangular sublattices of graphene, labeled A 
and B. The scattering at each interface conserves energy 
E = Elc and the component p y = hq of the momentum 
parallel to the interface, where the latter plays the role 
of channel index n (note, however, that due to spin and 
valleys, each value of q has degeneracy g — 4). In the left 
and right contacts, where Ul,r — >■ —00, a mode propa- 
gating towards the right (like the incident and transmit- 
ted ones) is proportional to (1,1), while a propagating 
mode moving towards the left (the reflected one) is pro- 
portional to (1, — 1) T . In the pump region, the scattering 
wave function with transverse wave vector q is given by 
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where A = ± is positive for electron-like and negative 
for hole-like quasiparticles, z ± = (±k x + iq)/y/k% + g 2 , 
and d=fc x = i\/[( e — U i)/L] 2 — q 2 is the electron's longi- 
tudinal momentum along the transport direction (+ for 
a quasiparticle moving towards the right, — for a quasi- 
particle moving towards the left). 

To facilitate the calculation, we first solve the scatter- 
ing problem of a single potential barrier, u\ for example, 
and then find the result for the double-barrier problem 
by composition of scattering matrices. The wave match- 
ing condition of continuity for a single barrier at x = 
is 
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Note that the phases from the infinite longitudinal wave 
vector in the contact regions (where Ul,r — >• 00) can be 
absorbed in the amplitudes T\ and t\. A similar set of 
equations can be written for a particle incoming from the 
right contact and transmitted to the left. The resulting 
scattering matrix has elements 
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where D = 



Xik x Li cos(k x Li) — ikFLism(k x Li), and 
2 The index Ai is +1 when 

1 



k x = ^{e-UxY/Ll -q 

the Fermi energy is above the onsite potential, or 
when it is below. The scattering matrix for the second 
barrier has the same form, with subscript 1 changed to 
2. The scattering matrix for the complete double barrier 
can then be calculated using the composition rule 
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We introduce the resulting S 2B into Eq. (lb), where 
(12 = ^1,2, and then take the weak-driving limit of Eq. 
pi. The index Ai = A2 then only depends on the position 
of the Fermi energy relative to the working point, and the 
longitudinal momentum in the pumping regions takes the 
form k x = \/k 2 F — <? 2 , where kp = e/L = E /(EpL) is the 
Fermi wave vector. Collecting all results for the graphene 
pump with heavily doped contacts (and L\ = L 2 = L), 
Eq. ^ yields 
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sin 2 (fcajL) [sin(2fcajL) — 2k x L cos(2k x L)] 
[(k x L)i + (qL)ism 2 (2k x L)] 2 

where the ± sign denotes whether the pump is doped 
with electrons (+) or holes (— ). This result not only 
applies to propagating modes (real momentum k x , \q\ < 
fci?), but also to evanescent modes, \q\ > kp, for which 
k x is imaginary. 

In contrast, a weakly driven graphene pump with un- 
doped leads has no incoming lead modes that become 
evanescent in the pump. This is because the Fermi mo- 
mentum in the pump at the working point is identical to 
the Fermi momentum in the contacts (Ul,r = 0). Taking 
the appropriate limit of the wave-matching results, each 
propagating mode then contributes a pumped charge 

_ _ k F L 2(qL) 2 cos(k x L) sin 3 (k x L) 



where \q\ < kp such that k x is real; there is no contribu- 
tion by modes with \q\ > kp. 

For normal pumps, the wave matching procedure is 
modified due to the different dispersion relations and ab- 
sence of the pseudospin degree of freedom; otherwise, the 
formalism remains unchanged. In terms of the longitudi- 
nal momentum K x of electrons in the leads, the pumped 
charge in each channel then takes the form 
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2k x L(k 2 x - K 2 X ) cos(2k x L) + (k 2 x + Kl) sin(2k x L) 

[4k 2 K 2 cos 2 (2k x L)^(k 2 ^K 2 ) 2 sm 2 (2k x L)} 2 

Both the heavily doped contact (Q^~°°) and undoped 
contact (Qq~°) normal cases follow from the above ex- 
pression by taking K x to infinity or k x: respectively. 

Equipped with Eqs. ([8|, (J9J), and the two appropriate 
limits of Eq. ( 10 ), we now can compare the results for the 
different settings. In all four cases, the pumped charge 
has a prefactor A u kp, indicating that pumping is propor- 
tional to the dimensionless driving strength A u and the 
pump's number N p = gkpW ji\ of propagating modes at 
the Fermi energy. The latter number includes a degener- 
acy factor g = 4 for graphene (accounting for two valleys 
and two physical spin states), and g = 2 for normal con- 
ductors (spin). By factoring out these two quantities, we 
obtain the dimensionless pumping response 
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which depends only on the system's scattering character- 
istics at a given energy. 

For short and wide systems (W ^> L), the dimension- 
less pumping response develops a quasi-continuous de- 
pendence on the transverse momentum q. In the two- 
dimensional contour plots of Fig. [2] we examine this de- 
pendence for varying carrier concentration (parameter- 
ized by the Fermi momentum kp), comparing the results 
for the two graphene setups (panels a and b) to those of 
the two normal pumps (panels c and d). For graphene, 
we only show the results for a Fermi energy above the 
working point; when the carriers in the central pump re- 
gion are changed from electrons to holes, the pumped 
current reverts sign if the leads are heavily doped [cf. Eq. 
d8|)], but remains the same if the leads are undoped [in 
any case the pumped current always changes sign when 
one reverts the pumping cycle]. 

In each panel, the dashed line |g| = kp delineates 
the border between propagating modes (|g| < kp) and 
evanescent modes (\q\ > kp). Pumps with undoped con- 
tacts cannot access such modes. For highly doped con- 
tacts, evanescent modes penetrate into the pumping re- 
gion, but only in the case of graphene [Fig. [21a)] their 
contribution is sizeable. This is especially true close to 
the Dirac point, where the evanescent modes dominate. 

This effect can be understood by considering the spe- 
cific conditions for electronic confinement in graphene, 




FIG. 2. (Color online) Momentum distribution of pumped 
charge per mode \q — Q q /(eA u N p ) as a function of mode in- 
dex q for varying carrier concentration (parameterized by the 
Fermi momentum /cf)- Blue and red represent opposite direc- 
tions of pumping (left to right or right to left). In graphene, 
the propagating mode with q — (normal incidence) cannot 
be pumped due to the Klein paradox. In the case of graphene 
with heavily doped leads significant pumping is possible due 
to the contribution of the evanescent modes (\q\ > /cf, delin- 
eated by the dashed line), which dominate around the Dirac 
point (/cf = 0). The other pumps can only drive current 
through the propagating modes. 



which are directly linked to the chirality of the charge 
carriers. Chirality conservation at the contact enables 
evanescent electrons to populate the graphene pumping 
region for modes within a window of width Aq ~ 1/2L 
around q = [18-20 . These evanescent modes con- 
tribute to pumping because they are sensitive to the on- 
site potentials Ui and have a finite amplitude at both 
contacts, so that charge transfer between them is pos- 
sible over a pumping cycle. It is noteworthy that, in 
contrast, the propagating mode with q = (normal in- 
cidence) never contributes to the pumping in graphene 
[both for doped as well as for undoped contacts, Fig. 
uib)]. This is a direct consequence of Klein tunneling 



-the transmission T n: 



1 is perfect at all energies, 



the mode is therefore insensitive to driving and cannot 
be pumped. 

Carriers in normal conductors do not display chirality, 
which entails that they can be easily confined. In partic- 
ular, for a normal conductor with heavily doped contacts, 
the large Fermi velocity mismatch suppresses the trans- 
parency of the contacts for all modes except those close to 
resonance with well-resolved energy levels. In Fig. [2|c), 
these levels are seen as narrow regions of finite pump- 
ing, along with a threshold kpL = ir/2 below which no 
pumping occurs. The pumping is directed, meaning that 
for a given orientation of the driving cycle, the pumped 
current has the same sign for all energies. The different 
confinement properties also entail that the contribution 
of evanescent modes to pumping is negligible at all ener- 
gies. 




or the precise position of the Fermi energy (as long as 
it is within less than El of the pump's Dirac point). It 
then acquires the universal dimensionless value 
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0.0288. (13) 



Due to its association to evanescent modes, this value can 
be interpreted as the pumping analogue to the minimal 
conductivity in stationary transport [18-20 . All other 
pumps have a vanishing pumping response at kp = 0. 



FIG. 3. (Color online) Pumped charge Q per cycle (normal- 
ized to the pumping strength A u and the number of propa- 
gating modes N p in short and wide pumps, W ^> L), as a 
function of the carrier concentration (parameterized by the 
Fermi momentum /cf). The plotted quantity is the pump- 
ing response \ u > introduced in Eq. (12), which is evaluated 
using Eq. ( |8|) ( graphene with highly doped leads, thick solid 
curve), Eq. (|9} (graphene with undoped leads, thick dashed 



curve), and the appropriate limits of Eq. (10) (normal pump 
with highly doped leads, thin solid curve, and with undoped 
leads, thin dashed curve). The inset focusses onto the re- 
gion close to charge neutrality {UfL <C 1), where evanescent 
modes in the graphene pump with highly doped leads allow 
for a finite charge transfer approaching the universal value of 
X u = 0.0288. 



The normal pumps become open when they are at- 
tached to undoped leads, and consequently in this case 
there is no energy threshold for pumping [Fig.[2rd)]. The 
sign of the pumped current is energy dependent, which is 
a generic feature of open pumps (including the graphene 
pump with doped leads). However, the contribution of 
evanescent modes vanishes identically since all incoming 
modes remain propagating in the pumping region. 

The total pumped charge can be characterized by sum- 
ming the mode-resolved result over all incoming modes, 



5E4 = irEa 



Q 



eA u N p 



(12) 



For short and wide systems (W ^> L), the sum can be ap- 
proximated by an integral over the continuous transverse 
momentum, J^ — >> (W/2n) J dq. 

The result for the four types of pumps is shown in 
Fig. [3j For large energies (kp > 1/L) the pumping re- 
sponse rises to 1/2 and 1/4 in the cases of undoped and 
highly doped leads, respectively. As a consequence of the 
contact-induced resonant tunneling subbands, the nor- 
mal pump with highly doped leads only operates above a 
finite carrier-concentration threshold. As highlighted in 
the inset, the evanescent electron pumping in graphene 
with highly doped leads results in a finite pumped charge 
at nominally vanishing charge-carrier density (kp = 0). 
In the considered limit W ^> L, the response is inde- 
pendent of the detailed system characteristics, like the 
system dimensions, aspect ratio (as long as it is large), 



IV. GATE VOLTAGE DRIVING 

In a realistic experimental pumping setup, the princi- 
pal driving parameters are not the onsite energies Ui, but 
gate voltages Vi (see Fig. fl]) which control the locally in- 
duced charge densities pi . The onsite energy is related to 
the charge density via the density of states, which differs 
between normal systems and graphene. In particular, the 
compressibility of graphene vanishes at the Dirac point 
[25] . In the following we will see that this enhances the 
response of a graphene pump. 

Because of screening, the translation between V and U 
in general requires a self-consistent treatment since the 
equilibrium charge density profile, and hence U(x), can 
be a complicated non-local function of V(x). The charge 
in the pump, which interacts with the whole gate struc- 
ture, equilibrates to the density p(x) that minimizes the 
total electrostatic energy. As a consequence, the onsite 
potential profile created by spatially piecewise constant 
gate voltages will acquire deviations from the piecewise 
constant model we have assumed for U(x) in the previ- 
ous sections. A second non-local screening contribution 
comes from the influence of the lead electrons on the den- 
sity close to the contact regions. This 'lead-doping' ef- 
fect is particularly relevant in graphene, which has more 
transparent contacts than a normal pump. As a result, 
the equilibrium electron density profile varies smoothly 
across the transparent contact from the lead's high den- 
sity to the lower density in the pump, which results in an 
inhomogeneous charge density (and hence onsite poten- 
tial) close to the contact. The equilibrium distribution 
depends on screening and contact details and was studied 
within an ab-initio approach in Ref. [26]. The screening 
of U(x) is expected to modify the results of the abrupt 
barrier model only quantitatively, since perfect normal 
transmission is preserved, although the angular profile 
away from normal incidence is narrowed. 

In view of these complications, we rely on the large ca- 
pacitance of the metallic gates to ignore the detailed local 
effects of screening in the pumping region, and account 
for non-local screening effects by relating the total charge 
rti under gate i = 1,2 to voltages Vj through a non- 
diagonal capacitance matrix Cij, rii(t) = J2jCijVj(t) 
(zero voltage is identified to charge neutrality). We then 
are in a position to assume that the charge density pi 



under each gate electrode is constant. The local onsite 
energy is related to the charge density by the integral of 
the density of states from the local position of the neu- 
trality point to the Fermi energy. Taking into account 
the different densities of states of graphene and normal 
conductors, one can then express the total charge ni in 
terms of the dimensionless onsite potential Ui = UiJEl 
and Fermi energy e = E /El used in our previous anal- 
ysis. For graphene rii(i) = eW(e — Ui(t)) 2 /(irLi), whilst 
for the normal case rii(t) = eW(e — Ui(i))/(27rLi). We 
then have 



whilst for a normal pump 



. W e 

L 2n 52j C iJ V o 
W e 



graphene 
normal 



(14) 



If we drive the pump with a weak adiabatic gate 
voltage cycle V-(t) = V^(0) + SVi(t), with V£(0) such 
that the working point is uniform, ui(0) = 1*2(0) = 
[and therefore ni(0) = 77,2(0) = ^(0)], we can apply 
the theory of the previous sections with a correspond- 
ing cycle in the onsite potential parameter space ui{t) = 
(dui/dVj)SVj(t). The resulting areas enclosed in the u 
and V parameter spaces are then related by the Jaco- 
bian A u = Ay det (dui/dVj). Specifically 
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for the normal pump. Note the n(0) in the denominator 
of A gr . This implies that the onsite energy response to ex- 
ternal gate driving in graphene diverges at the neutrality 
point, and can be seen as a direct consequence of the van- 
ishing electronic compressibility, n oc (du/dV) -1 —> 0. 
Furthermore, we observe that the pumped current in the 
weak driving regime A u <C 1 is suppressed by the cross 
capacitance between the electrodes (since they diminish 
det Cij). The optimal driving is still achieved by maxi- 
mizing the area Ay , which for the typical case of sinu- 
soidal driving corresponds to a phase shift of tt/2 between 
the two gates. 

The pumped current at driving frequency v is given 
by I = vQ = evA u N p x u (kFL). We can express / as a 
function of Ay and the initial charge under each gate, 
n(0) = Y, j C ij V j (0), by using 
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This gives for graphene 
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When driving via gate voltages at fixed strength Ay, the 
pumping efficiency W j Ay diverges in graphene close to 
its incompressible neutrality point (n(0) — >• 0), while it 
vanishes in the normal case. Recall, however, that by 
assumption SVi(t) <C "^(0), so that Ay det C^ <C n 2 (0). 
Therefore, although the pumping efficiency for weak driv- 
ing may diverge in graphene, the current itself will not. 



V. PUMPING LIMITATIONS 

As a general rule, the pumped current increases with 
increasing pumping frequency v for small v. In the adi- 
abatic regime studied here it does so linearly in u : but 
this optimal scaling typically becomes sublinear as one 
approaches the non-adiabatic regime, hv ~ El, [27 pos- 
sibly with superimposed oscillations depending on the 
model [28] . For a small graphene pu mp of L = 0.1/im, 
this adiabatic frequency ceiling is around ^o.i/im ~ 1.6 
THz, which drops down to ^i Mm ~ 0.16 THz for a larger 
L = 1/im pump. This is a very high frequency if we com- 
pare it to the one of normal pumps, ^o.i/im = 14 GHz 



and v\. 



0.14GHz. 



Regarding the maximal magnitude of adiabatically 
pumped current, both normal (fabricated e.g. on 
GaAs/AlGaAs heterojunctions) and graphene pumps are 
comparable, providing an estimated 10 to 500 nA for typ- 
ical setups with many electrons in the pump. Assuming 
W/L = 6, C - 1 aF and V(0) - 10 V (which corresponds 
to an n(0) of around 62 electrons under each gate) and 
a driving of SV ~ IV, we have a maximum current of 
around 225 nA (725 nA) for graphene (normal) pumps 
of length L = 0.1/im, which drops down to 23 nA (7 
nA) for L = 1/im. The advantage of graphene in this re- 
spect becomes most noticeable in the few-electron regime 
n(0) ~ 1. As an example, the maximum current in an 
L = 0.1/im graphene pump at V(0) ~ 10 SV ~ 1 V 
is still in the 0.8 nA range due to the evanescent-mode 
contribution, while the normal pump is effectively inop- 
erative in this range, due to the first subband threshold. 
A further advantage of graphene versus normal pumps is 
that the ballistic transport regime faborable for pumping 
which is assumed throughout this work is much more easy 
to achieve (especially at high temperatures) in graphene 
than in semiconducting heterostructures. 

We finally discuss thermal fluctuations, which gener- 
ally degrade quantum pumping (possible exceptions are 
Refs. |221 EDI )• At the very least, if T < hv/k B , 
one would just observe an effective thermal smearing of 
x(kpL) in Fig. pi since in this regime thermal fluctua- 
tions can be considered static within each driving cycle 
[3T] , Above this temperature, thermal fluctuations will 
quickly suppress quantum pumping, since quantum co- 
herence within a single pumping cycle is destroyed. Since 



v has a ceiling given by adiabaticity, this allows us to es- 
timate the maximum operating temperature of adiabatic 
pumps by Tl — E^/kB-, which is around 76 K for a 
small L = 0.1/im graphene pump, as opposed to 7.6 K 
for a comparable normal pump. This thermal advantage 
arises as a result of the density of states. It is also one 
of the reasons for the temperature robustness of other 
transport effects in graphene. 



VI. CONCLUDING REMARKS 

In summary, we have investigated fundamental and 
practical aspects of adiabatic quantum pumping which 
distinguish graphene-based systems from equivalent se- 
tups involving conventional two-dimensional electron 
gases. Because of the unique properties of graphene (in 
particular, the chirality of the charge carriers), evanes- 
cent modes can contribute significantly to the pump- 
ing, especially when the system is operated close to the 
charge-neutrality point. For the case of short and wide 
pumps, the evanescent pumping regime is characterized 
by a universal value of the dimensionless pumping re- 
sponse. This value does not depend on the width or 



length of the pump. It is also largely independent of 
temperature, as long as T < E^jhs^ since the evanes- 
cent mode pumping response is quite flat within energies 
~ El of the Dirac point. In normal pumps evanescent 
modes only give a negligible contribution. 

In practical terms, the vanishing electronic compress- 
ibility of graphene at the Dirac point enhances the re- 
sponse of graphene-based pumps to driving via external 
gate potentials. This presents a clear performance advan- 
tage for graphene pumps if they are operated in the few 
electron regime. Furthermore, graphene-based quantum 
pumping promises an enhanced robustness against ther- 
mal effects, as already known from stationary transport. 
For the same reason, the regime of adiabatic driving ex- 
tends to higher frequencies than in normal pumps. As 
for other graphene-based electronic applications, these 
attractive features are further enhanced by the long co- 
herence time and high mobility of charge carriers in 
graphene. 
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